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The title material has a quasi-one-dimensional electronic structure and is of considerable interest 
because it has a metallic phase with properties different from a simple Fermi liquid, a poorly 
understood "insulating" phase, and a superconducting phase which may involve spin triplet Cooper 
pairs. Using the Slater-Koster approach and comparison with published band structure calculations 
we present the simplest possible tight-binding model for the electronic band structure near the 
Fermi energy. This describes a set of ladders with weak (and frustrated) inter-ladder hopping. In 
the corresponding lattice model the system is actually close to one-quarter filling (i.e., one electron 
per pair of sites) rather than half-filling, as has often been claimed. We consider the simplest 
possible effective Hamiltonian that may capture the subtle competition between unconventional 
superconducting, charge ordered, and non-Fermi liquid metal phases. We argue that this is an 
extended Hubbard model with long-range Coulomb interactions. Estimates of the relevant values of 
the parameters in the Hamiltonian are given. NMR relaxation rate experiments should be performed 
to clarify the role of charge fluctuations in Lio.9MoeOi7 associated with the proximity to a Coulomb 
driven charge ordering transition. 

PACS numbers: 71.27.+a, 74.40.Kb, 71.30.+h, 71.45. Lr 

I. INTRODUCTION 

The electronic properties of quasi-one-dimensional materials (Q1D) are of particularly interest because of the 
possibility that they may exhibit properties characteristic of Luttinger liquid rather the conventional Fermi liquid 
metallic state seen in three dimensional metals. These systems may display Luttinger liquid (LL) behavior above a 
temperature at which the thermal energy fcgT is larger than the interchain hopping energy t±— . In a LL charge and 
spin are carried by independent collective excitations instead of the spin- 1/2 charged quasiparticles present in a Fermi 
liquid. Photoemission experiments can probe the existence of a LL through the observation of low energy spinon 
and holon branches and a power-law suppression of the density-of-states at the Fermi energy. Possible realizations 
of a LL have been investigated extensively in recent years. Examples are the Q1D cuprate materials, SrCu02 (Ref. 
H) and Sr2Cu03 (Ref. y), quasi-lD organic crystals^, carbon nanotubes^^, GaAs channels^ and more recently 
one-dimensional Au chains deposited on Gc(001) surfaces^. 

Apart from the intrinsic interest in finding a LL in actual materials, quasi-onc-dimensional materials often exhibit 
low temperature broken symmetry states including superconducting, charge and/or spin density waves and Peierls 
phases. These instabilities of the metallic phase occur at a crossover temperature scale at which the system effectively 
becomes three-dimensional and so long range order can occur at low but finite temperatures. The physics associated 
with this crossover which goes from the pure one-dimensional to the three-dimensional system as temperature is 
reduced and understanding the possible competing ground states at low temperatures remains a formidable theoretical 
challenge. 

Another example of a material which may exhibit one-dimensional physics is the blue bronze Lio.9MoOeOi7 which 
shows some behavior consistent with a LlJ^ at sufficiently high temperatures. Angular Resolved Photoemission 
Spectra (ARPES) shows that quasiparticles do not exist in the system since a power law suppression in the DOS is 
found at the Fermi energy and dispersing spinon and holon branches are seen in the spectral density. At sufficiently 
low temperatures there is a transition to an insulating-like or semiconducting state from which superconductivity 
occurs at even lower temperatures. The low temperature "insulating" state is poorly understood. The experimental 
data pose challenging questions to address apart from the LL behavior observed above the crossover temperature 
scale. In order to understand the mechanism of superconductivity, it is important to determine the symmetry of the 
Cooper pairs. It is also important to characterize the "insulating" phase seen below T m ~ 25 K. Finally, it is necessary 
to characterize the nature of the excitations in the metallic state, in particular, if the material is a "bad metal" with 
incoherent excitations or not. 

Here we introduce a realistic microscopic model on a lattice which can capture the essential physics observed. This 
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FIG. 1: (Color online) A speculative temperature-pressure phase diagram relevant to Lio.gMo60i7. It assumes that at ambient 
pressure the system is close to a quantum critical point for an insulating charge ordered phase. Fluctuations associated with 
this phase cause the resistivity to increase with decreasing temperature below T m . The pressure dependence of this curve is 
taken from Fig. 3 (a) in Ref. [ll|. Associated with the quantum critical point (QCP) there is a crossover from a quantum 
critical metal (or "bad" metal) to a Fermi liquid metal at a temperature scale To- The associated coherence temperature goes 
to zero as the QCP is approached. Similar ph ysic s has been observed in recent studies of a two-dimensional extended Hubbard 



model consists of weakly coupled ladders describing the d-orbitals in the Mo atoms with strong on-site and off-site 
Coulomb repulsion. Based on this model we establish that Lio.gMoOgOi? should be regarded as a nearly quarter-filled 
system and not as a nearly half-filled system. This means that in the absence of longer range Coulomb interactions the 
systems would always be a metal. Thus the long range part of the Coulomb repulsion should not be neglected and can 
lead to charge ordering phenomena and a Wigner-Mott insulator. The model contains rich physics since it contains 
low dimensionality, strong Coulomb interactions, geometrical frustration, and charge order frustration. We discuss 
possible ground states of the proposed model in different limits. From the model proposed, experimental observations 
and theoretical considerations we speculate that Lio.gMo060i7 is close to a quantum critical point dominated by 
charge ordering fluctuations. A possible schematic phase diagram is shown in Fig. [1] 

The present paper is organized as follows. In Sec. |TT] we briefly revise key experimental observations in 
Lio.gMo060i7. In Sec. IHII we discuss the one-electron tight-binding parameters appropriate for Lio.gMo060i7 based 
on a Slater-Koster parameterization and comparison with DFT calculations of the band structure. Sec. I IVI introduces 
the Coulomb repulsion energies which are relevant for the description of the title material. This leads to a minimal 
strongly correlated model which consists on weakly coupled ladders comprising the full three dimensional crystal 
structure. Finally, in Sec. [V]wc discuss the theoretical implications of the model obtained and possible consequences. 



We now review some key observed properties of Li 9 Mo 6 Oi7, with a particular emphasis on deviations from the 
behavior seen in conventional metals and charge ordered insulators. 



Recent measurements of the resistivity along the three crystal axes yield o, resistivity anisotropy of: p a : p^ : p c — 
80:1:1600 at T = 300 K and p a : p b : p c ~ 150: 1: 1600 at T = 4.2 K^. Much smaller anisotropies have been 
observed by other author a 16 ' 17 . However, care must be taken in measuring the conductivity parallel to the chains 
because it can be reduced significantly if the current path is not strictly parallel to the chains*^ These are consistent 
with observed anisotropies in the superconducting upper critical fields^. Hence, these measurements corroborate the 
quasi-one-dimensionality of the system in its electronic properties. 




II. BRIEF REVIEW OF EXPERIMENTS 



A. Electronic anisotropies 
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B. The bad metal 

Transport properties of a wide range of strongly correlated electron materials exhibit certain features that are 
distinctly different from the electronic properties of elemental metals^. These unusual properties arise from emergence 
of a low energy scale which defines a temperature scale T (often in the range 10-100 K) above which quasi-particles 
do not exist and the material is referred to as a bad metal. Signatures of this crossover from a Fermi liquid at low 
temperature to a bad metal at temperatures above T include: (i) the resistivity, Hall coefficient, and thermopower 
have a non-monotonic temperature dependence, (ii) with increasing temperature the intralayer resistivity smoothly 
increases to values much larger than the Mott-Ioffe-Regel limit (h 2 a/e ~ 1 mficm) corresponding to a mean-free 
path comparable to the lattice constant ai^, (iii) at temperatures of order To the thermopower has values as large 
as ks/e ~ 80//V/K, and (iv) above To the Drude peak in the frequency dependent conductivity collapses and the 
associated spectral weight shifts to higher frequencies. 

It is important to realise that if quasi-particles do not exist then the one-electron spectral function does not have 
dispersive features, i.e., it is completely incoherent. This means that above To the notion of a band structure and a 
Fermi surface have no meaning. All of the above features in the transport are captured by a dynamical mean-field 
theory (DMFT) treatment of lattice Hamiltonians such as the Hubbard model and Anderson lattice models. It is 
generally believed that a small To is associated with proximity to a Mott insulator or to a quantum critical point^ 2 -. 
An alternative picture of bad metals can given in terms of hard core bosons 2 ^. 

Lio.gMogOi7 exhibits many of the above signatures of bad metallic behavior. 

Resistivity. The magnitude of the resistivity parallel to the chains is about 100 /ificm at 100 K . 14 ' 15 This is 
comparable to that observed in other "bad metals" such as optimally doped cuprate superconductors! 19 ' 21 For a 
quasi-onc-dimensional Fermi liquid the conductivity can be written as — Ae 2 £/(ac)h where £ is the mean- free path 
along the chains and a and c are the unit cell dimensions perpendicular to the chain direction. This leads to an 
estimate of t ~ 106 at 100 K. 

Thermoelectric power. The thermopower is positive in sign and a non monotonic function of temperature, increasing 
from a value of approximately 20 /iV/K at low temperatures to a maximum of about 30 /uV/K at 50 K and then 
decreasing slowly to a value of about 10 //V/K around room temperatures^ 2 -. Recent measurements^ 3 - found the 
magnitude of the thermopower exhibited a significant sample dependence, with values as large as 200 /iV/K, and 
was interpreted in terms of bipolar transport in an almost perfectly compensated metal. All of the above behavior 
contrasts to what is observed in a simple Fermi liquid in which the thermopower is approximately linear in temperature 
for temperatures less than Tp 2 ^. Furthermore, in contrast to experiments the band structure suggests that the 
thermopower should be negative since the system has an electron-like Fermi surface, with two bands close to one- 
quarter filling. 

Frequency dependent conductivity. The optical conductivity at both 10 K and 300 K is relatively flat as a function 
of frequencies of the scale of tens of meV— , suggesting the absence of a Drude peak. The total spectral weight at low 
frequencies is quite small, corresponding to about one charge carrier per 10 unit cells if one integrates up to about 
0.1 cV. For electric field polarisations perpendicular to the chain direction there is a clear energy gap of about 0.4 
eV, at both 10 K and 300 K. Hence, the material has the very unusual property that it appears to be conducting for 
transport parallel to the chains and insulating for transport perpendicular to the chains. 

Hall effect. The transverse conductivity o~ a b increases by a factor of about 60 as the temperature decreases from 300 
K to 25 K.— The Hall coefficient increases by several orders of magnitude from (1 — 3) x 10 _9 m 3 /C to about 10 _6 m 3 /C 
as temperature decreases from 200 K to 2 Ki 16 ' 23 This is inconsistent with a simple Fermi liquid metal in which the 
Hall coefficient should be independent of temperature. The value obtained for Lio^MoOgOiy assuming a Drude 
model, Rh = —l/n\e\ is Rh ~ —2 x 10 -9 m 3 /C, assuming n = 1.9 electron charge carriers per unit cell. This value is 
consistent with the absolute experimental value at T = 250 K although with the opposite sign.— The T-dependence 
of Rh shares some similarities with the observed behavior in the organic charge transfer salt (TMTTF^AsFg. In this 
material, the Hall coefficient increases by a factor of about 100, as the temperature decreases from 200 K to 100 K, 
at which there is a charge ordering transitional The Hall coefficient has been calculated for a set of weakly coupled 
Luttingcr liquids with umklapp scattering associated with one-half filling 2 ^. Similar results are expected near one- 
quarter filling. There are small power law temperature-dependent corrections to the high temperature non-interacting 
value. The bandwidth sets the scale for this temperature dependence. In summary, the Hall effect is inconsistent with 
both a quasi- ID Fermi liquid and a Luttingcr liquid picture even at the qualitative level. 

Thermal conductivity. The Lorenz ratio is about 10-30 times the Fermi liquid value, implying a gross violation 
of the Wiedemann-Franz law^ This has been interpreted in terms of a Luttingcr liquid picture with spin-charge 
separation. The Hall Lorenz ratio increases from about 100 to 10 5 times the Fermi liquid value as the temperature 
decreases from 300 K to 25 

Angle- Resolved Photoemission Spectroscopy (ARPES). This measures the one-electron spectral density and 
result s 10 ' 29 have been interpreted in terms of Luttinger liquid and non-Fermi liquid pictures. The ARPES data 
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displays characteristic holon and spinon branches as well as the characteristic suppression of the single particle den- 
sity of states (DOS) of a one-dimensional system signalling the destruction of quasiparticles and of Fermi liquid 
behavior. However, the exponent a characterizing the DOS suppression depends on T varying between: a w 0.6 for 
T < 200K and a « 0.9 at T = 300K in contrast with the T-independent a in a one band Luttinger liquid and has 
been attributed to renormalization of a due to interaction of charge neutral critical modes associated with the two 
bands crossing the Fermi energy. The scaling of the spectral density with temperature: A(k,Lu) = T v A(vk/T,u>/T), 
violates the one-band LL scaling relation: r\ — (a — 1), since experimentally r\ = 0.56 instead of rj < 0. This unconven- 
tional quantum critical scaling indicates the presence of quantum fluctuations which can mask the pure LL behavior. 
Interestingly no warping of the Fermi surface has been observed in ARPES yet. 

Nernst effect. Recent measurements have been interpreted in terms of bipolar transport in an almost perfectly 
compensated metal with close to equal numbers of electrons and holes^ The magnitude of the Nernst signal is four 
orders of magnitude larger than the value given by Behnia's simple Fermi liquid expression^ 

Magnetoresistance. This exhibits a number of unusual properties, for all temperature a 14 ! 16 . For magnetic fields 
parallel to the chains there is a small negative magnetoresistance, for all current directions, suggesting suppression of 
the "insulating" state. For magnetic fields and currents both perpendicular to the layers (the c-axis direction) there 
is a huge positive magnetoresistance. At fields of 10 Tclsa it can be as large as 50 and 500-fold at temperatures of 
50 K and 3 K, respectively— This is non-classical as there is no Lorentz force since the magnetic field and electric 
current are parallel. Similarily, unusual behaviour has also been seen in a wide range of other strongly correlated 
low-dimensional metals. 



C. The "insulating" state 

Some properties suggest the possibility of an insulating or semiconducting-like phase at low temperatures. However, 
it should be stressed that one does not see activated behaviour (i.e. clear evidence for an energy gap). It is observed 
that the resistivity is a decreasing function of temperature from the superconducting transition temperature T c ~ 1 
K up to about 20 K. The resistivity then increases approximately linear with increasing temperature up to room 
temperature. The simplest possible explanation would be that there is a metal- insulator transition around 20 K. 
However, this is inconsistent with several experiments we discuss below. First, the minimum in the resistivity occurs 
at T m at different temperatures for different current directions, ranging from 15K < T m i n < 30A". Increasing the 
pressure to 20 kbar T m decreases from about 30 to 10 K (Figure 1), and the magnitude of the resistivity decreases 
significantly-ii. Hence, it is not clear that the "insulating" state exists at high pressures. We also note that high 
magnetic fields parallel to the chain direction reduce the low temperature resistivity, which can be interpreted as a 
destruction of the "insulating" stated. 

The fact that the resistivity is a decreasing function of temperature above the superconducting transition is rather 
unusual and puzzling since one normally sees a direct transition from a metallic phase to a superconducting phase. 
However, there are other cases involving quarter-filled organic charge transfer salts where a superconducting state 
occurs close to a charge ordered insulator [see the Table in Ref. Hi] and the inset of Fig. 2 in Ref. [32[. Other 
quasi-one dimensional materials exhibit a resistivity with a similar temperature and pressure dependence similar to 
that summarised in Figure 1. For example, Per2M(mnt)2 [M=Pt,Au] is a CDW insulator at ambient pressure but 
above 0.5 GPa the resistivity has a non-monotonic temperature dependence^ 

An important question is whether the "insulating" state is a Charge Density Wave (CDW) driven by a Fermi surface 
instability with a partially gapped Fermi surface. Such CDWs are observed in other quasi-one-dimensional materials 
such as K0.3M0O3 and the transition-metal dicalchogenide compound 2H-NbSe2^. The structural instabilities as- 
sociated with the CDW can be clearly seen in x-ray scattering. However, high resolution x-ray scattering, neutron 
scattering^, and thermal expansion^ experiments on Lio.gMogOi7 observe no structural instability such as one would 
expect to be associated with charge ordering phenomena. However, observing a structural instability driven by Fermi 
surface nesting requires a sufficiently large electron-lattice coupling which may not be present. Other Mo compounds 
such as the blue bronze K0.3M0O3 which is quasi-two-dimensional show jumps in resistivity typical of a more con- 
ventional CDW accompanied by a Peierls transition. A possible interpretation of the observations in Lio.gMogOiy is 
that an unconventional electronically-driven CDW occurs which is not detectable with structural analysis data. For 
example, the "insulating" state may be a D-Density Wave (DDW) state, which has nodes in the energy gap. Such a 
state has been proposed as the low-temperature state of a-(BEDT-TTF)2KHg(SCN)4 (Ref. l37h and of the pseudogap 
state in the cuprates^. Since in a DDW state there is no modulation of the charge density in real space it is hard to 
detect the associated symmetry breaking or long range order. 

It is important to determine the possible nature of the magnetic interactions if they exist in the "insulating" 
phase. The temperature dependence of the magnetic susceptibility has a Curie contribution from a small number 
of magnetic impurities (about 10~ 4 per conduction electron)^ Subtracting the Curie contribution, the remainder 
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is weakly temperature dependent, between about 10 K and 200 K, characeristic of the Pauli paramagnetism of a 
Fermi liquid metal. Hence, there is no sign of the energy gap that is normally seen in CDWs and charge ordered 
insulators. Pauli susceptibility data^ lead to: x(0) ~ 2.8 x 10~ 6 consistent with previous data^ 9 .. At low temperatures 
the specific heat capacity has a term that is approximately linear in temperature with a coefficient 7 ~ 6 mJ / (mol 
K 2 )i 39 i 40 . More recent data find that the specific heat coefficient at temperatures right above the superconducting 
transition temperature is: 7 ~ 1.6 mJ/ (mol K 2 ). From the Fermi velocity obtained from band structure calculations 4 ^ 
hvp = 3.7eVA, one finds a bare density of states at the Fermi energy due to the two Mo(d xy ) one-dimensional bands: 
D(cf) = 1-9 states/(eV cell). This unrenormalized bare DOS is three times smaller than the DOS obtained from the 
experimental^ specific heat slope: 7 ~ 1.6 mJ /mol K 2 ). Experimental data lead to a Sommcrfcld- Wilson ratio^: 
R = 47r 2 fc|jx(0)/(3((7^B) 2 7) ps 2, indicating substantial electronic correlation effects^. 

The fact that the charge transport properties (resistivity and Hall coefficient) show "insulating" behaviour with 
decreasing temperature while the thermodynamic properties (specific heat and magnetic susceptibility) show Fermi 
liquid properties is puzzling. One possible explanation is the latter are associated with low energy spin excitations 
and not charge excitations. 

D. Superconducting state 

The transition temperature T c can vary between about 1-2 K depending on sample purity. Substituting Li ions with 
K and Na ions led to a reduction in T c and a correlation between T c with the residual resistivity^ Such a sensitivity 
is characteristic of unconventional superconductor, i.e., non s-wave pairing4^ Increasing the pressure to 20 kbar T c 
increases from about 1.8 to 2.5 K.— The upper critical field for magnetic fields parallel to the chains may be above 
the Chandrasckhar-Clogston-Pauli paramagnetic limit for a spin singlet superconductor— suggesting the possiblity of 
a spin triplet state. 

E. Isoelectronic materials 

It might be expected that the materials Ao.gMoeOi7 with the alkali metals A=K,Na,Tl have similar properties 
and much can be learned from comparisons with A=Li. However, it turns out that these materials have a slightly 
different crystal structure, with separated metal-oxygen layers leading to a significantly different electronic structured 
Specifically, they turn out to have three partially filled d-block bands and a quasi-two-dimensional band structure 
and Fermi surface. They undergo CDW instabilities due to a hidden Fermi surface nesting^ For A=K the CDW 
transition occurs at 120 K, associated with a structural transition increasing the unit cell dimensions four fold, as seen 
by X-ray scattering and electron diffraction.— At the transition about 50 per cent of the charge carriers are gapped 
out. For A=Na the CDW transition occurs at 80 K, and opening of the energy gap on two of the bands crossing the 
Fermi energy has been seen in angle- resolved photoemission spectroscopy (ARPES)4£ 

III. TIGHT BINDING BAND STRUCTURE 

In order to introduce the simplest strongly correlated model to describe the electronic properties of Lio.gMo060i7, 
we first analyze the band structure published earlier using the extended Hiickel method^ and the Local Density 
Approximation (LDA) version of Density Functional Theory (DFT)^ 1 -. 

These calculations allow extracting the nominal valence of the compound, LiMoeOi7. The band structure shows the 
existence of two d xy bands crossing the Fermi energy and two d xz /d yz bands which arc filled with two electrons leading 
to six electrons in the unit cell. The d xz and d yz bands are shifted away from the d xy bands and the Fermi energy 
becoming complctly filled due to hybridization with the neighboring O atoms. This effectively leads to two isolated d xy 
bands crossing the Fermi energy which contain two electrons. Since there are four Mo atoms per unit cell, this implies 
that each Mo chain atom has 1.5 electrons with one electron in the d xz /d yz orbitals and half an electron in the d xy 
orbitals. This corresponds to a chemical valence for LiMogOi7 of Li +1 (Mo 4,5+ )2(Mo' 6+ )4(0 2 ~)i7 where Mo' denote 
atoms which are not in the chains. On the other hand, the doped compound: Li +0 - 9 (Mo 4 ' 55+ )2(Mo' 6+ )4(O 2_ )i7, has 
the Mo orbital with supressed electron density: Mo(4d 145 ) so that there is less than 0.5 electrons (about 0.45) per 
d xv orbital of Mo. 
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TABLE I: Position vectors of the four Mo atoms*- inside the unit cell of Lio.gMo60i7 in units of the monoclinic a, b and c 
unit cell vectors (a — 12.762A, b — 5.523A and c = 9.499A). The angle between the a and c axes of the monoclinic crystal is 
P = 90.61° 



Moi 




Ra/o 




Moi 


0.9939 


0.25 


0.23356 


Mo 4 


0.16635 


0.25 


0.9206 


Mo'i 


0.00613 


0.75 


0.7664 


Mo' 4 


0.8337 


0.75 


0.07938 



A. Slater-Koster parameters 

We parameterize the model using the Slater-Koster tight-binding approac h 49 ' 50 . We use 2p x ,2p y ,2p z orbitals as 
the minimum basis set for describing O and Ad xy orbitals for describing Mo. Using distances obtained from the X-ray 
crystal structure 5 ** we analyze the hopping amplitudes between Mo and O atoms and the direct hopping amplitudes 
between the Mo atoms. In a second stage, we single out the important hopping amplitudes needed for an effective 
tight-binding model Hamiltonian describing the d xy orbitals of Mo. 

The distances and displacement vectors needed to define the cosines between the various atoms are shown in Table U 
and [XT] together with estimated hopping amplitudes obtained from the Slater-Koster approach. The hoppings between 
the Mo and O atoms are given by: 

tp^ v = m(V3l 2 V pd a + (1 - 2l 2 )V pdv ) (1) 
t Pv d xv = l(V3m 2 V pda + (1 - 2m 2 )V pdn ) 
tpxd xy = lmn(V3V pd a - 2V P cl-k) 



where I, m, n are the direction cosines of the vector from the p orbital to the d orbital and the parameters: 

29.5 cV 



u MoO 



(2) 



Vpd7r ~ y/ 2 



13.6 eV 

~~7i 

a. 



"MoO 

where the distances are in A and the parameters are appropriate for Mo. According to the above equations one finds 
that only the t p d m hopping amplitude between Mo and O atoms in a chain is non-zero and the corresponding value 
shown in Table [ill Also the hopping amplitude between different chains in a ladder through an O atom is zero with 
the Slater-Koster approach. The direct hopping between two d xy orbitals maybe estimated from: 

td^y = 3l 2 m 2 V dda + {I 2 + m 2 - M 2 m 2 )V Mnx (3) 

with the corresponding parameters: 

213.3 eV 

Vdda = --Jg (4) 

a MoMo 



Vddl! — 



115.2 eV 
~d? 

M0M0 



where again the distances are all in A. The nearest neighbor intra-chain hopping between Mo is then tM 0l Mo' = Vdd-w 
at a distance d = 3.725A and between Mo atoms at d = 3.69A in different chains: ijwoiMo 4 = VddS = [50[. The 
hopping amplitude between Mo atoms on different ladders is non-negligible: tMo 1 Mo' 1 = Vdd-K = 0.018eV. 

The largest contribution to the interaction between neighboring Mo atoms in the chain comes from hopping through 
intermediate O atoms. A straightforward estimate of this hopping amplitude from perturbation theory gives: 

t 2 

t(Mo - Mo) = Pydxy = 0.9 eV, (5) 
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TABLE II: Distances between molybdenum and oxygen atoms*i and hopping amplitudes in Lio.9MoeOi7. Hoppings are 
estimated based on the Slater-Koster approach. The most relevant nearest-neighbor Mo-Mo and Mo-0 distances within the 
main ladder and between different ladders are displayed. Hopping amplitudes are between 4d xy orbitals Mo atoms and between 
4:d xy and 2p y orbitals for Mo-0 within a chain. Distances are given in A and hopping in eV. The labeling of the atoms follows 
Ref. M (see Fig. 2). 

Atom-atom Distance(A) Hopping (eV) 



Moi 


-Mo' 4 (in-chain) 3.725 


0.1606 


Moi- 


-On (in-chain) 1.873 


1.515 


Moi 


-M04 (inter-chain) 3.6756 


0. 


Moi 


-Oi (inter-chain) 1.873 


0. 


Moi- 


-Mo'i (inter-ladder) 5.7655 


0.018 



which leads to a half-bandwidth of 1.8 eV for the chain, a value which is large compared to the LDA-DFT value 
of about 0.9 — 1 eV (Ref. l4ll ). Such a discrepancy may be related to the breakdown of the perturbative form 
since:|eMo — eo| ~ tp y d xy using the values CMo(4d) = —11.56 eV and eo(2p) = —14.13 eV. We will discuss in more 
detail below how the effective nearest-neighbor hopping between Mo atoms needs to be reduced by a factor of two 
with respect to the Slater-Koster hoppings in order to correctly reproduce the LDA-DFT bandwidths. 

From the Slater-Koster analysis above we would arrive at the conclusion that the zig-zag Mo chains are then 
decoupled and only a non-zero inter-ladder hopping exists. However, DFT band structure calculations do indicate the 
existence of inter-chain hopping which is, however, much weaker than the effective in-chain hopping. This is because 
the direct Moi-Mo'4 bonding in the chain is essentially of the 7r-type whereas the interchain M01-M04 bond is of the 
weak £-type since the chains are on top of each other. Therefore, the inter-chain coupling should come from the small 
5-type coupling between the d xy orbitals of the Mo atoms. Ladders are coupled through a direct hopping connecting 
Mo d xy orbitals. The interladder Moi-Mo'i bonding is of the 7r-type. On the other hand, hopping processes via two 
intermediate O atoms should be very small. This is because one of these intermediate oxygens is an apical O, (O4 in 
Fig. 2 of Ref. l5ll ) weakly hybridized to Moi via the O4 and O's p-orbitals. It is important to note that interchain 
and interladder hoppings should be non-zero in order to recover the dispersion along the c-direction found in DFT 
caculations. A three parameter Slater-Koster approach^ to the bands would lead to an interchain hopping amplitude: 
tMoxMoi = VddS — Vdda/6 = —0.046 cV which has the opposite sign to the 7r-typc of bond between different ladders 
^MoiAfo'j = 0.018 eV. We will see below how the opposite signs between these two hopping amplitudes are important 
in a tight-binding model to capture the LDA-DFT band structure close to the Fermi energy. 



B. Effective tight-binding model for Lio.oMosOir 

We have performed a tight-binding calculation in which we include the four d xy orbitals of the Mo atoms and 
integrate out the O atoms which are assumed to enter indirectly via the Mo — Mo hopping amplitudes. We keep 
the most relevant hopping amplitudes needed for the description of the band structure. The simple tight-binding 
model for the d xy orbitals is compared to previous LDA-DFT band structure calculations and the hopping amplitudes 
extracted. 

The position of the atoms in the unit cell are expressed as: Ri = Ri a a.+ Ribh + Ri C c and the momentum wavevector: 
k = k a a* + fc(,b* + k c c* , all refered to the unit cell coordinate system. A translation vector of the lattice reads: 
Rn = n a a + nt,b + n c c and the Bloch functions: 

|* Moi (k)) = -L^VMR^+Rmo, )|* MOi (r-R n -R Moi )) (6) 

^ R„ 

For the effective in-chain interaction we have: 

(* Moi (k)|lf|* Mol (k)) = (* Mo4 (k)|ff|vI/ Mo ,(k)) = -2tA(k a ,k c )cos(^-), (7) 

where A(k a ,k c ) = e i27r ^ fca ' fljf °4 RM °i^ +k ^ R Mo> 1 r mo-l)c) ^ interaction between chains in the same ladder is de- 
scribed through matrix elements as: 

(y Moi (k)\H\y Moi (k)) = _t ±e *-(a-c) e *k.(R Mo4 -R Moi ) = - t± j*S± (g) 
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(0,7i/b) 



(0,0) 



(7C/C.0) 



FIG. 2: (Color online) Band structure for the lower two bands obtained from the diagonalization of the 4x4 tight-binding 
matrix defined in the text using t = 0.5 eV, t± — —0.024 eV, t' = 0.036 eV. These parameter values are chosen to produce 
band dispersions comparable to LDA-DFT calculations^. The band dispersions along the b and c directions are plotted. The 
horizontal line denotes the Fermi energy. 



where S± = 0.17a — 0.31c, and similarly for the Mo^-Mo^ interaction: 

(tf A/oi (k)|Jf|tf A/oi (k)) =-f ± e" ik '^ = (3Wk)|if|3Wk)}* 
Two ladders in two neighboring unit cells are coupled through the matrix element: 

(* A/oi (k)|#|* Moi (k)) = _t' e *V k ' (RM °i~ RMoi) _ t ' e * (a-b) e *.(R Moi -R Moi ) 



-2iV k51 cos( — ) 
v 2 ' 



(9) 



(10) 



where <5i = a — b/2 + R-Moi — R-Moi = 0.01a + 0.53c. On the other hand the M04 — Mo' 4 hopping amplitude between 
ladders is weaker since they are at a larger distance and is neglected: 



(*M 04 (k)|If|*M<(k)) =0. 
The final tight-binding 4x4 Hamiltonian to be diagonalized reads: 



(11) 



/ 



H(k) 







-t±e 



-2i'e ik<51 cos^) 



-t 1 e 



-2t'e- lk Sl cos( k ^) -2tA(k a , k c ) cos(^) 



-2tA(k a , k c ) cos( 




2 > 



-2tA(k a ,k c )cos(^) \ 




V -2tA(k a ,k c )cos(^) 







-t±e 



In the absence of intcrladder hopping (t 1 = 0), the four bands: 

e(k) = ±t± ± 2tcos( 



kb, 



— tj_e 



— zk-<5 j 







J 



(12) 



which recovers the expected dispersions of the uncoupled two-leg ladders, as it should. 

Diagonalizing the Hamiltonian we obtain the band structure along the b and c-directions shown in Fig. [2] The 
effective hopping between the nearest-neighbor Moi-Mo'4 atoms in a chain is taken to be t — 0.5 cV in order to 
reproduce the DFT bandwidth of the two lowest Mo d xy energy bands. We use the hopping amplitudes t' = 0.036 cV 
a factor of two larger than the one obtained from Slater-Kostcr in Table [TT] and we take the hopping between chains 
in the ladder opposite to t 1 and of magnitude: t± = —0.024 eV consistent with the type of bonding (tt vs. S bonding). 
This set of parameters captures the correct magnitude of the band dispersions along 6, the dispersion of the two lowest 
bands in the c-direction due to the combined coupling between ladders, t' , and the interchain hopping amplitude, t±. 
The magnitude of the values extracted from the present analysis t± and t' of about 0.03 eV are consistent with recent 
preliminar estimates^. 

The Fermi surface obtained from the diagonalization of the 4x4 Hamiltonian is shown in Fig. [3J This tight-binding 
Fermi surface is close to the DFT Fermi surface (cf. Fig. 6. in Ref. [4ll ) except for the latter the two Fermi surfaces 
touch at the zone boundary. The filling is n — 1.9, with the warping of the most filled band opposite to the less filled 
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FIG. 3: (Color online) Fermi surfaces obtained from the diagonalization of the 4x4 matrix using t = 0.5 eV, t± = —0.024 eV, 
t' = 0.036 eV and filling n = 1.9. The Fermi surfaces corresponding to the two separate bands are denoted in blue and red. 
Note that the two sheets have opposite warping. 




q(l/b) q(l/b) q(l/b) 

FIG. 4: (Color online) Wavevector dependence of the three non-interacting charge susceptibilities. The three boxes correspond 
to decreasing temperature from left to right. The three susceptibilities are associated with the two low energy bands crossing 
the Fermi surface and shown in Fig. [3] with t = 0.5 eV, t± = -0.024 eV, t' = 0.036 eV and filling n = 1.9. The interband 
susceptibility X12 displays an instability at low temperatures due to significant nesting between these two Fermi sheets. 



band. We note that the opposite signs of t' and t± found in the Slater-Koster approach is essential for capturing the 
opposite warping of the two Fermi sheets found in DFT. The Fermi surface obtained from the diagonalization of the 
4x4 Hamiltonian is shown in Fig. [3] 

Fermi surface nesting instabilities. The tendency of the system towards nesting instabilities can be analyzed by 
computing the bare charge susceptibility, 

Xw(q) = em(k + q) _ em , (k) • ( 13 ) 

between the different bands, e m (k). From a simple inspection of the Fermi surface in Fig. [3] one observes that since 
the two Fermi sheets corresponding to the two bands are oppositely warped the nesting condition is optimally satisfied 
with a nesting vector which is about q w (0,7r/6, 0) between the two different bands. Evaluating Xmm'(q) shown in 
Fig. |H a divergence is observed in the interband susceptibility: Xi2(q) when decreasing the temperature but not in 
the intraband contributions. 



C. Quarter filling or half filling? 

From a band structure point of view the apparent half filling arises because of the zig-zag structure of the chains 
which leads to a folding of the Brillouin zone. However, the key point (at least in terms of strong electronic correlations) 
is that for LiMoOeOi7 there would be one electron per two Mo ions. In this sense the material should be viewed 
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as being close to one-quarter filling. This is important because it means that long-range Coulomb interactions are 
required to produce an insulating state. We also note that both ARPES and LDA give kp ~ ir/2b and this only 
corresponds to half filling when b is the nearest-neighbour distance. 



D. Anisotropy in the conductivity 

In a simple Fermi liquid picture the ratio of the conductivity parallel and perpendicular to the chain directions is 
independent of temperature and given by 



<?b_ ^ ft 
<Ja ~ V t' 



(14) 



The estimates in Table J II T| suggest an anisotropy ratio of about 100 which is consistent with that observed experi- 
mentally in References 

MM, and[H but much larger than reported in References [l?] and [l6l . 



IV. EFFECTIVE STRONGLY CORRELATED HAMILTONIAN 



An effective Hamiltonian which can capture the electronic properties of Lio.gMo060i7 is now presented. We argue 
that since the system is close to quarter-filling not only the onsite U but also the off-site Coulomb repulsion V should 
be included to account for charge ordering phenomena. This minimal model is an extended Hubbard model which 
reads: 

H = H a + Hu, (15) 
where Hq is the non-interacting tight-binding Hamiltonian which reads: 

H = ~t E (cSUW + cc) -t ± J2 ic$L£ff + c.c.) 

- -t'Vfc (,)t c {l+1) +cc)-t ,s r(c m c {l+1) +cc) 

(16) 

where the first term describes the kinetic energy of a single ladder and the second term the hopping processes between 
ladders. The index I denotes one of the two chains in a ladder, a labels a specific ladder and i runs from 1 to N, the 
number of sites in the chains. The parameters t, t±, denote hopping amplitudes within a ladder whereas t' denotes 
the hopping amplitude connecting nearest-neighbor ladders. One only needs three hopping amplitude parameters for 
the description of the band structure. The Coulomb interactions are encoded in Hjj which reads: 



Hu = U £ n«na + E fn<MO m +^E "S^f + V E^Sffi + n% +1 n^) + V" E n%n% +2 

cx,l,i<T a,/,z a,/,z a,/,z ct,/,i 

+ ^E«xi +1 +^"E«S«- 

(17) 

The parameters U, V, V±, V and V" denote intraladder Coulomb interactions, whereas V'", V"" and V"" denote 
Coulomb interactions between nearest-neighbor ladders. A schematic representation of the parameters entering the 
model proposed are displayed in Fig. [S] In principle, the model keeps the essential intra-ladder and inter-ladder 
Coulomb repulsion energies that may be relevant for Lio.gMoOgOiy. Estimates of the values of Coulomb repulsion 
energies entering the Hamiltonian are provided and discussed below. 



A. Model parameter values 



In Table Mil we show the hopping amplitudes and Coulomb parameters estimated for model (fT5")) . From the 
comparison to the DFT-LDA calculations we found that the nearest-neighbor hopping is t = 0.5 eV, the hopping 
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FIG. 5: (Color online) Schematic representation of the model Hamiltonian (|15[) to describe L10.9M06O17. The hopping and 
Coulomb energies entering the model describing electrons in the d xy orbitals of Mo are displayed. A schematic view perpen- 
dicular to the b — c plane is shown noting that the planes of the ladders in the real crystal are tilted with respect to the b — c 
plane. The rectangle corresponds to the unit cell in the b-c plane. 



between chains in the same ladder is t±_ = —0.024 cV and between nearest-neighbor chains in different ladders: 
t' = 0.036 eV. We note that for the case in which t' = — 3tj_/2, we have that the two bands cross at k c = in 
agreement with the DFT calculations. When the same sign is used t' = 3t±/2, then the two bands do not cross at 
k c = 0. So it is essential that t±_ and t' do have opposite signs to capture the appropriate warping and dispersions. 
The on-site Coulomb repulsion is about U = 6.4 eV estimated from constrained-DFT— calculations. 

The unscreened Coulomb interaction between electrons in nearest neighbor Mo atoms within a chain is estimated 
to be V = 2 cV, which is comparable to the bandwidth and therefore relevant unless complete screening inside the 
crystal occurs. These values are large enough to lead to charge ordered ground states since the critical value V c in a 
quarter-filled chain is about: V « It. Screened values of V arc reduced to about V — 0.2 eV, assuming a dielectric 
constant of ew 10eo due to the screening of the rest of the crystal^ 

In a similar quarter- filled ladder— compound NaV2 0s charge ordering phenomena due to the off-site Coulomb 
repulsion has been found. From DFT-based calculations the parameters are estimated to be: t\\ ~ 0.17 eV, t± ~ 3t, 
and U ~ 2.8 eV. The distance between vanadium ions along the chain is approximately 3.8A and the nearest-neighbor 
Coulomb repulsion V has been estimated^ to be V « It consistent with the zig-zag type charge ordering^ found in 
these systems. 

Experimental observations combined with the present discussion suggest that in Lio.gMogOi7, it is desirable that 
the Coulomb repulsion V ~ 0.5 — 1 eV so that the system is placed near the QCP to a charge order transition driven 
by V (cf. Fig. [1]). This requires that the long-range Coulomb repulsion is only partially screened and is important in 
order to understand the electronic properties of Lio.gMogOi7. The relevance of incomplete Coulomb screening to the 
electronic properties of quasi-one-dimensional systems has already been pointed out in the case of TTF-TCNQ 55 i 56 
crystals. 

We now discuss the relevance of further neighbor Coulomb interactions. The Coulomb repulsion between chains in 
the same ladder should be similar to the nearest-neighbor V, since the distance between two neighboring Mo atoms in 
the same chain is 3.725 A and is 3.675 A between atoms in different chains, so V± ~ V. This parameter is important 
for stabilizing a zig-zag type of charge ordered state in the ladder. Comparing the distances between different Mo 
atoms one realises that further neighbor interactions are relevant to the model. In fact, the Coulomb repulsion between 
electrons on Mo atoms on different zig-zag chains of the same ladder (see Fig. [5]), V' is comparable to V and V± since 
interladder distances are between 5.256 and 5.366 A. Since the next-nearest neighbors distance between two Mo is 
5.52 A due to the zig-zag nature of the chains, we have that V" is not be negligible. On the other hand the nearest 
distance between Mo atoms in nearest ladders is 5.77 A so the Coulomb repulsion between neighboring ladders V'" is 
comparabe to V". The off-site Coulomb repulsion included in the model is cut-off at values of the Coulomb repulsion 
at which: V"" ~ V/2 ~ t, assuming V ~ 2t. Coulomb energies between farther distant neighbors are neglected since 
they are smaller than the bare kinetic energy of the electrons. 

Hence, a minimal model should contain the nearest-neighbor Coulomb interactions and possibly longer range 
Coulomb interactions between electrons in further distant Mo atoms as represented in Fig. [5] Accurate screening 
calculations in L10.9M06O17 should be performed in order to pin down the values of the off-site Coulomb repulsion and 
settle the relevance of the longer range Coulomb repulsion in L10.9M06O17. Below we discuss the minimal strongly 
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TABLE III: Estimates of the Coulomb and hopping parameters parameters for Lio.9MoeOi7. The on-site Coulomb repulsion, 
U, is obtained from constrained-DFT calculations^. We have assumed: V = It appropriate for having the material in the 
proximity to a charge order QCP (see Fig. [!}. We have assumed a 1/d decay of the further off-diagonal Coulomb parameters 
with d, the distance between electrons in different sites. The Coulomb repulsion energies and hopping amplitudes in this table 
can be identified in the schematic structure of Fig. \E\ All energies are given in eV. 

t__t_L t' U V V x V V" V" V"" 
0.5 -0.024 0.036 6.4 1 1 0.7 0.66 0.65 0.53 

correlated model (fT5|) in some specific limits, for which there are known results. 

The estimated hoppings and Coulomb parameters of the model are summarized in Table Hill 

V. UNDERSTANDING THE MODEL HAMILTONIAN 

In principle, the model (|15p is effectively two dimensional consisting of zig-zag ladders which are weakly coupled 
through t' and there is also Coulomb repulsion between electrons in different ladders. The relevance of the inter ladder 
hopping interactions is settled by the temperature scale at which actual experiments are undertaken. As a first step 
in the understanding of the full complexity of the model we discuss the physics of isolated chains and ladders for 
which much more theoretical work is available. Therefore, we first consider the model (|15|) in different limits assuming 
independent ladders (t' = and V'" = V"" = 0) so we are left only with t, t± and the Coulomb repulsion energies: 
U, V, Vx and V. 

A. What is the ground state of the ladder model at 1/4 filling and U — > oo? 

We consider quarter-filled ladders assuming that doubly occupied sites are forbidden (U — > oo). This model has 
been studied numerically using DMRG^ and analytically through weak coupling RG and bosonization techniques^ 
in different parameter ranges which we now discuss. 

1. tx = Vx = V = 

The two chains comprising the ladder are completely decoupled and the model maps onto the t — V model of 
spinless fermions at half filling. The ground state is a Luttinger liquid for V < 2t and an insulator with long range 
charge order for V > 2tM- Even in the presence of the long range Coulomb interaction decaying as 1/d, where d is 
the distance between electrons in different sites, it is found that the Wigner state with one electron at every other 
site of the lattice is the ground state at one-quarter filling^ 5 - when t — s- 0. This is satisfied in this case indicating 
that charge ordered states are non-frustrated by the long range part of the Coulomb interaction. Also recent DMRG 
calculations on two-leg ladders with tx = t show an insulating phase with a charge gap which is interpreted as a result 
of dimerization of the rungs with one electron localized on each rung of the ladder—. 

2. tx = V = 0, V ± / 

This ladder model has been studied by bosonisation and RG approaches^. There are two charge modes, total 
and difference. For < V < 2t, the difference mode can be gapless. The other is gapped in the presence of a small 
Vx- This is an homogeneous insulating phase consistent with DMRG results performed for tx/t < 1.—. When both 
V = Vx > 2t the insulating phase has zig-zag charge order with wave vector (ir,ir). The role of Vx is crucial in 
locking the charge order waves running along the two independent when these are decoupled. Interestingly, even at 
1^ = and for any finite Vx the system displays a charge gap which has an exponential dependence with Vx as in the 
half-filled Hubbard model leading to an homogeneous insulator. 
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3. V ^ 

This interaction frustrates charge order and should produce a metallic state when V' ^ V ^ V±, even when 
V In this limit, the model can be mapped onto a classical Ising model with frustrated interactions which leads 

to a disordered ground state. A quartcr-hllcd one-dimensional extended models with nearest, V, and next-nearest 
Coulomb interactions, V' has been analyzed using DMRG which shows a metallic state when V' ~ V/2. 

B. Triplet superconductivity 

Quantum Monte Carlo calculations on a nearly quarter-filled model of weakly coupled chains with onsite Coulomb 
repulsion interaction U only and small U — 2t display /-wave (spin triplet) superconducting tendencies due to 2k /- 
CDW instabilities. The /-wave symmetry is related to the fact that the electronic modulation is of about four 
lattice spacings: Q ps 2k f w tt/2 close to quarter- filling^. Such behavior has been confirmed by a random phase 
approximation (RPA) analysis in weakly coupled quarter-filled chains with Coulomb interactions up to third nearest 
neighbors in the presence of interchain Coulomb repulsion^. It is found that for U = 1.7t, t± = 0.2t and under 
the condition: V' + V±_ ~ U/2, triplet /-wave superconductivity wins over the d-wave channel. This is because the 
presence of V± enhances the 2fc/-CDW instabilities of the isolated chains. Triplet superconductivity pairing has also 
been encountered in weak coupling RG calculations on weakly coupled chains in the presence of both intrachain and 
interchain Coulomb repulsion^ for moderate values of the interchain Coulomb repulsion. 

For the two-dimensional quarter-filled extended Hubbard model it is found that introducing frustrating charge 
interactions on the square lattice (e.g., along one diagonal) destabilises the stripe charge ordered insulating phases 
producing a charge ordered 3-fold symmetric metallic stated at large V and at finite-!/. In the limit of U >> t 
a metallic 'pinball' liquid state in which no doubly occupied sites ('pins') can occu r 12 i 67 . Melting the 3-fold state 
produces an 'f-wave' spin triplet superconducting state. How this particular symmetry emerges can be seen from 
direct inspection of the 3-fold charge ordering pattern in real space which consists of placing electrons on the closest 
sites not coupled via the Coulomb interaction to a given occupied site, avoiding the off-site Coulomb repulsion. Thus, 
the /-wave pairing found in the isotropic triangular lattice is analogous to the d xy -~w&ve pairing found on the square 
lattice^! and results from electrons avoiding the strongest nearest neighbor Coulomb repulsion. 

C. Doping the ladder away from one-quarter filling 

Since Lio.gMogOi7 is slightly doped away from one quarter-filling it is worth considering this situation. Assuming 
that t± — V = and Vj_ ^ as above there is a metallic state if the chemical potential is of the order of the gap of 
the quarter-filled system^ based on weak couping RG. This is because a commensurate- incommensurate transition 
occurs since Umklapp processes are suppressed as: 4fc/ ^ |G| where G is a reciprocal vector of the lattice. On the 
other hand, superconducting fluctuations in the doped quarter-filled ladder are found to be of the singlet d-wave type 
but dominated by the Akf-CDW correlations. However, including longer range Coulomb interactions can change this 
picture and induce /-wave triplet pairing. It remains an open question to understand how the presence of V can 
influence the superconducting tendencies of the system in ladders doped away from quarter-filling. 

VI. CONCLUSIONS 

The title quasi-one-dimensional material displays an intriguing competition between insulating, superconducting 
and "bad" metallic behavior. Besides the mechanism of superconductivity, the nature of the "insulating" phase and 
the unconventional metallic properties are poorly understood. No evidence of a structural transition has been found 
accompanying the occurrence of the "insulating" phase. Under pressure the insulating phase is suppressed giving way 
to conventional metallic behavior below a low temperature crossover scale. 

In order to understand these phenomena, we have derived a minimal strongly correlated model for determining 
the low energy electronic properties of Lio.gMogOiy. The one-electron part of the Hamiltonian is obtained based 
on a Slater-Koster approach and comparison with DFT band structure calculations. The tight-binding Hamiltonian 
consists on three hopping parameters only: t, t± and t' which capture the dispersion of the two Mo(d xy ) bands 
along the 6-direction crossing the Fermi energy, the opposite warping of the two Fermi surface sections, and the weak 
dispersion of the two bands in the c-direction. These are the main features found in full band structure calculations. 
The real space tight-binding Hamiltonian describes zig-zag ladders weakly coupled by the small inter-ladder hopping, 
t'. In the lattice model we note that the system is close to quarter-filling although the bands are nearly half- filled. 
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This is a result of the band folding associated with the four atoms in the unit cell arising from the zig-zag structure 
of the ladders. 

A reinterpretation of the physics of Lio.gMogOi7 as a nearly quarter-filled system instead of a nearly half-filled 
system is suggested from the model introduced. Both the onsite U and long range Coulomb repulsion are found to 
be relevant since estimated values can be comparable to the bandwidth of the material. Based on these estimates 
and experiments under pressure we suggest that Lio.9MoeOi7 is close to a charge ordering transition driven by the 
Coulomb repulsion. Based on the sensitivity of the material to external pressure we argue that many of the anomalies 
observed may arise due to the proximity to a Quantum Critical Point (QCP). The "bad" metal behavior is attributed 
to quantum criticality: a metal with unconventional excitations arising from the charge fluctuations occurring at all 
length scales around the QCP. Associated with the QCP is the existence of a crossover temperature T* below which 
coherent excitations and Fermi liquid behavior occurs on the metallic side of the transition. This scale is suppressed 
as the QCP is approached. Quantum fluctuations associated with charge order occur in the proximity to the QCP 
which can lead to the enhancement in the resistivity below T m leading to a resistivity minimum which has been found 
in other quasi-one-dimensional materials close to CDW instabilities, such as Per2M(mnt)2 [M=Pt,Au].— 

Spin triplet superconductivity can arise in systems in which charge fluctuations associated with a nearby charge 
ordered phase dominate. In general, in ladder systems explored close to quarter-filling, zig-zag charge order correla- 
tions are strongly enhanced even by moderate interchain Coulomb repulsion, V±, as previous works have shown. In 
Lio.gMogOi7, the fact that V± = V and the presence of next-nearest neighbors Coulomb interaction and longer range 
interactions along the chains due to the zig-zag structure can favor spin triplet superconductivity. 

Experiments probing the existence of large charge fluctuations near the QCP are desirable in order to understand 
the enhancement of the resistivity with decreasing temperature. X-ray diffraction find no structural changes associated 
with the upturn of the resistivity at T m suggesting that a purely electronic mechanism plays a major role. From the 
T-dependence of the 1/TiT-NMR spectra one could find whether anomalous broadening or splitting of spectral finest 
occurs when lowering the temperature below T m in analogy to TMTTF2ASF6 (Ref. lo9h , in which the charge ordering 
transition detected by NMR is not accompanied by a change of the structure. This kind of experiments could clarify 
the role played by charge ordering phenomena driven by the long range Coulomb repulsion as proposed here. 
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